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Abstract.  Tf.is  paper  proposes  an  improved  high-resolution  technique  that  incorporates  the  forward  (original)  array  together 
with  its  complex  conjugated  backward  version  to  achieve  superior  performance.  This  is  in  turn  realized  by  averaging  the 
forward  covariance  matrix  and  the  backward  one,  and  using  it  in  conjunction  with  the  MLJSIC  technique.  This  paper  also 
presents  the  asymptotic  analysis  in  terms  of  first-order  approximations  of  the  mean  and  variance  of  the  null  spectrum  estimator. 
It  is  shown  here  that  in  an  uncorrelated  scene,  the  bias  corresponding  to  this  scheme  is  exactly  half  of  that  associated  with 
the  conventional  one  while  maintaining  the  same  variance.  Furthermore,  the  bias  expression  in  the  new  scheme  is  used  to 
obtain  a  resolution  threshold  for  two  uncorrelated,  equipowered  plane  waves  in  white  noise  and  the  result  is  compared  with 
that  jbtained  by  the  conventional  scheme  for  the  same  source  scene. 

Zusammenfassung.  Dieser  Beitrag  schlagt  ein  verbessertes.  hochauflosendes  .Spektralanalyseverfahren  vor,  bei  dem  das 
Vorwarts-Array  und  das  zugehorige  konjugiert-komplexe  Riickwarts-Array  kombiniert  werden,  urn  besseres  Verhalten  zu 
erzielen.  Dies  wird  dadurch  erreicht,  da|5  der  Mittelwert  von  Vor-  und  Ruckwartskovarianzmatrizen  gebildet  und  anschliePend 
der  MUSIC-Algorithmus  verwendet  wird.  Auperdem  wird  eine  asymptotische  Analyse  von  Mittelwert  und  Varianz  des 
Nullraum-Spektralschatzers  vorgestellt.  Es  wird  gezeigt,  dap  bei  unkorrelierter  Umgebung  Mittelwert  und  Vaiianz  genau  unt 
den  Faktor  zwei  gegeniiber  dem  herkommlichen  Verfahren  reduziert  werden  konnen.  Die  Ausdriicke  fiir  den  Schatzfehler 
werden  dartiberhinaus  verwendet,  um  eine  Auflosungsgrenze  fiir  zwei  unkorrelierte,  im  weiBen  Rauschen  eingebetlete  ebene 
Wellen  gleicher  Leistung  anzugeben.  Dieses  Ergebnis  wird  mit  der  Grenze,  die  das  herkommliche  Verfahren  fiir  gleiche 
Eingangssignale  liefert,  verglichen. 

Rnume.  Cet  article  propose  une  technique  a  haute  resolution  amelioree  qui  incorpore  la  sequence  avant  ( originate )  avec  sa 
version  conjugue  complexe  arriere  pour  atteindre  des  performances  superieures.  Ceci  est  realise,  a  son  tour,  en  calculant  la 
moyenne  de  la  matrice  de  covariance  avant  et  arriere  et  en  I'utilisant  en  conjunction  a\cc  la  technique  MUSIC.  Cet  article 
presente  egalement  I'analyse  asymptotique  en  termes  des  approximations  d'ordre  un  de  la  moyenne  et  de  la  variance  de 
I'estimateur  spectral  sans  effet.  II  est  montre  ici,  que  dans  une  scene  sans  correlation,  le  biais  et  la  variance  correspondant 
a  cette  methode  sont  exactement  la  moitie  de  ceux  associes  avec  la  methode  conventionnelle,  En  plus,  les  expressions  de 
biais  de  la  nouvelle  methode  sont  utilisees  pour  obtenir  un  seuil  pour  deux  ondes  planes  non  correlees  et  a  energie  egale 
dans  du  bruit  blanc  et  le  resultat  est  compare  a  celui  obtenu  par  la  methode  conventionnelle  pour  la  meme  scene  de  source. 

Keywords.  Underwater  signal  processing,  detection  and  estimation  of  multiple  signals. 


1.  Introduction 

The  eigcnstructure-based  high-resolution  techniques  for  estimating  the  arrival  angles  of  multiple  plane 
wave  signals  have  been  of  great  interest  since  the  well-known  works  of  Pisarenko  [14],  Schmidt  [16]  and 
others  [2-4,  7,9,  11,  17],  These  methods,  in  general,  utilize  certain  eigenstructure  properties  resulting  from 
the  special  structure  of  the  sensor  array  output  covariance  matrix  for  planar  wavefronts  [16]  to  generate 
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spectral  peaks  (or  equivalently  spectral  nulls)  along  the  actual  directions  of  arrival  whenever  the  sources 
are  at  most  partially  correlated. 

When  the  exact  ensemble  average  of  the  array  output  covariances  are  used,  all  these  methods  result  in 
unbiased  values  (i.e.,  zero  for  the  null  spectrum )  along  the  true  arrival  angles,  irrespective  of  signal-to-noise 
ratios  (SNRs)  and  angular  separations  of  sources.  However,  when  these  covariances  are  estimated  from 
a  finite  number  of  independent  snapshots,  these  techniques  exhibit  deviations  from  their  ensemble  average 
values.  These  deviations  depend  on  the  specific  scheme  under  consideration  together  with  the  SNRs  and 
other  signal  and  array  specifications.  All  these  taken  together  determine  the  resolution  capacity  of  the 
technique  under  consideration. 

This  paper  proposes  an  improved  high-resolution  technique  without  compromising  on  the  size  of  the 
original  array.  In  addition  to  the  forward  (original)  array,  this  new  scheme  makes  use  of  a  complex 
conjugated  backward  version  of  the  original  one  to  achieve  superior  performance.  This  is  in  turn  realized 
by  averaging  the  forward  covariance  matrix  and  the  backward  one,  and  using  this  averaged  matrix  in 
conjunction  with  the  MUSIC  estimator.  A  detailed  performance  analysis  of  this  new  scheme,  when 
covaiiances  estimated  from  a  finite  sample  size  are  used  in  place  of  their  ensemble  averages,  is  also 
presented  along  with  first-order  approximations  of  the  mean  and  variance  of  the  null  spectrum  estimator. 

The  organization  of  this  paper  is  as  follows;  For  clarity  of  presentation,  the  conventional  scheme  is 
summarized  and  the  proposed  one  is  described  in  Section  2.  Using  results  derived  in  Appendix  A,  Section 
3  presents  the  first-order  approximations  of  the  mean  and  variance  of  the  null  spectrum  estimator 
corresponding  to  this  new  scheme  and  the  conventional  MUSIC  scheme.  The  bias  expressions  are  then 
used  to  obtain  a  resolution  threshold  for  two  uncorrelated,  equipowered  plane  wave  sources  in  white 
noise;  and  this  result  is  compared  to  the  resolution  threshold  in  the  conventional  case  [8]  for  the  same  scene. 


2.  Problem  formulation 

Let  a  uniform  array  consisting  of  M  sensors  receive  signals  from  K  narrowband  sources  M|(f), 

U;(t) . UK(t),  which  are  at  most  partially  correlated.  Furthermore,  the  respective  arrival  angles  are 

assumed  to  be  0, ,  0^, . , . ,  6^  with  respect  to  the  line  of  the  array.  Using  complex  envelope  representation, 
the  received  signal  .x,(t)  at  the  ith  sensor  can  be  expressed  as  [12] 

,v,(t)=  y  u^(f)e  (1) 

k  1 

Here  the  inter-element  distance  is  taken  to  be  half  the  wavelength  common  to  all  signals  and  n,{l) 
represents  the  additive  noise  at  the  ith  sensor.  It  is  assumed  that  the  signals  and  noises  are  stationary, 
zero-mean  circular  Gaussian'  independent  random  processes,  and  further,  the  noises  are  assumed  to  be 
independent  and  identical  between  themselves  with  common  variance  cr'.  Rewriting  ( 1 )  in  common  vector 


'  A  complex  random  vector  i  is  defined  to  be  circular  Gaussian  if  its  real  part  x  and  imaginary  part  >■  are  jointly  Gaussian  and 
their  joint  covariance  matrix  has  the  form  |5.  10] 


f: 


V 


where  :  =  x-rjy.  When  :  has  zero  mean,  its  covariance  matrix  is  given  by  =  f:  [( x  +  jyK.v'  -  jy'  I]  =  V  +  j  W.  Clearly,  E | =  0. 

Here  onwards  4'.  4'  stand  for  the  transpose  and  the  complex  conjugate  transpose  of  4,  respectively. 
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notation  and  with  coi^-  t!  cos  0^,  k  =  \  ,2 . K,  we  have 

_  K 

x'(l)-\M  ^  u^(t)a(w^)  +  n(/), 

K  I 

where  x^(t)  is  the  M  x  1  vector 

r'(n  =  [.x,(t),.v:(f) . -Vm(/)]^ 

and  a(wt)  is  the  normalized  direction  vector  associated  with  the  arrival  angle  i.e.. 

«(tuj  =  ^[l,e  '•“‘,e  '-“‘,...,e  '' ''  "“’‘f. 

V  M 

The  array  output  vector  j:'(t)  can  further  be  rewritten  as 


where 


.r'( =  Au{t) 


n(  f)  =  [u,(n,  u.(/) . U/c(t)]\  n(l)  =  [n,('K  «:( t) _ _  ristl  t)]^ 


A  =  v'M[a(w|),  a((0:) . aiw^  )].  (5b) 

Here  A  is  an  MxK  matrix  with  Vandermonde-structured  columns  {M>K)  of  rank  K.  From  our 
assumptions  it  follows  that  the  (forward)  array  output  covariance  matrix  has  the  form 


where 


E[/it)/\t)]  =  AR[A*  +  CT-I, 


R[^E[u(i)u\t)] 


represents  the  source  covariance  matrix  which  remains  as  nonsingular  so  long  as  the  sources  are  at  most 
partially  correlated.  In  that  case  AR[A*  is  also  of  rank  K  and  hence,  if  A,  A- >  ■  ■  •  ^  Av»  and 
^1 ,  ,  Pm  are  the  eigenvalues  and  the  corresponding  eigenvectors  of  /?'  respectively,  i.e., 

R'  =  1  \iP,Pl  (7) 

;  -  1 

then  the  above  rank  property  implies  that  A^  +  ,  =  A^  ,;  =  •  •  =  A^,  =  cr’  and  the  eigenvectors  corresponding 
to  these  equal  eigenvalues  are  orthogonal  to  the  direction  vectors  associated  with  the  true  arrivals,  i.e., 
P ,a(wi,)  =  0,  i  =  K  +  \,  K  +2, . . . ,  M,  k  =  \,2, . . . ,  K.  Hence  the  K  nulls  of  P(a>)  given  by 

Q(w)=  I  ||3;a(a;)p=l-  V  l/3;a(w)|’  (7a) 

k-’K  *\  )<  -  I 

correspond  to  the  actual  directions  of  arrival  [16]. 

To  improve  the  performance  of  the  above  estimator,  in  addition  to  the  forward  (original)  array,  we 
propose  to  make  use  of  the  complex  conjugated  backward  array  of  the  original  one.  Let  x*'{i)  denote  the 
complex  conjugate  of  the  backward  array,  i.e.. 


jr'’(/)  =  f,xt,(t),.vt,  ,(f) . .vf(f)]'  =  4H'’(t)+ii'’(t), 


v-i  r,  jv:  , 


I 

1»-/  k  I 


(8) 
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where 

and 

n'iO  =  [nit,(0,  . . . ,  (8b) 

It  follows  that  the  backward  array  output  covariance  matrix  has  the  form 

R^=E[x^{t)x^Ht)]  =  AR':.A"  +  a-I,  (9) 

where 

Rl  =  E[u^(t)u^\i)].  (9a) 

Averaging  the  forward  covariance  matrix  and  the  backward  one  together,  we  define  the  forward/back¬ 
ward  (f/b)  covariance  matrix  as 

R^k(R^+R^)  =  AR^A'  +  a^I,  (10) 

where 

+  (10a) 

It  is  easy  to  show  that  whenever  the  rank  of  J?^  is  at  least  K  -  I,  the  same  is  true  for  R^  and  moreover 
in  that  case  R^  is  of  rank  [12, 13].  Hence  the  eigenvalues  of  R  satisfy  -  >kK^.,  =  Kk*2  = 

•  •  ■  =  Am  Consequently,  as  in  (7a)  the  K  nulls  of  Qiu>)  given  by 

Q(oi)=  I  \pla(co)\~=l-  I  (11) 

(k  =  K  +  l 

will  correspond  to  the  actual  directions  of  arrival.  Here  p,,  -  are  the  eigenvectors  of  R 

corresponding  to  the  eigenvalues  A|,A2,...,Am;  i-C-, 

M 

R=l 

/=i 

Notice  that  the  above  lower  bound  on  the  rank  condition  (rank  Ri  =  K  -\)  physically  corresponds  to 
one  coherent  arrival  among  the  K  sources;  and  further  the  full  rank  property  for  the  equivalent  source 
covariance  matrix  R^  in  that  case  implies  that  the  coherent  source  has  been  essentially  decorrelated  by 
the  simultaneous  use  of  the  fc"  vard  and  backward  arri  H<. -^’ever,  the  main  advantages  of  this  modified 
scheme  do  not  become  apparent  until  its  performance  .  uation  is  completed  in  the  more  realistic 
“data-only  known”  case. 

So  far  we  have  assumed  that  an  ensemble  average  of  the  array  output  covariances  are  available.  Usually, 
these  exact  averages  are  unknown  and  they  are  estimated  from  the  array  output  data.  Often  this  is  carried 
out  for  the  unknowns  of  interest  using  the  maximum  likelihood  procedure.  For  zero-mean  M-variate 
(circular)  Gaussian  data  x\t„),  n  =  \,2, . . . ,  N  in  (5),  with  unknown  M  x  M  covariance  matrix  R\  the 
maximum  likelihood  (ML)  estimate  of  the  covariance  matrix  is  given,  referring  to  [1],  as 

i  xUt„)x^'(t„). 

N  „„t 

Using  the  invariant  property  of  the  maximum  likelihood  procedure,  ♦be  corresponding  estimates  S*’ 
and  S  for  the  unknown  matrices  R*'  and  R  can  be  constructed  from  by  the  same  rule  that  is  used  in 
constructing  R*'  and  R,  respectively,  from  R\  i.e., 

5  =  |(S'+.S"), 
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with 

N  „  =  i 

In  what  follows  we  study  the  statistical  properties  of  these  estimated  covariance  matrices  and  their 
associated  sample  estimators  for  direction  finding. 


3.  Performance  analysis 

3.1.  Main  results 

This  section  examines  the  statistical  behavior  of  the  proposed  forward/backward  scheme  and  derives 
expressions  for  the  bias  and  the  resolution  threshold  for  two  equipowered  uncorrelated  sources.  These 
results  are  substantiated  by  simulation  studies  and  comparisons  made  with  similar  results  obtained  in  the 
conventional  case  [8].  Toward  this  purpose,  consider  the  eigen-representation 

S  =  ELt  (12) 

for  the  ML  estimate  of  the  f/b  matrix  R,  where 

E  =  . L  =  diag^[r,, /, . Ik,  Ik* . ,hi],  and  EE*  =  1^, 

which  satisfies  i  =  1,  2, . . . ,  M  for  uniqueness.  Here  the  normalized  vectors  Ci,  Cj . Ck  are  the 

ML  estimates  of  the  eigenvectors  ^2,  -  •  ■ ,  0k  of  *  respectively  [5].  Similarly,  /,,  Ji . Ik  are  the 

ML  estimates  of  the  K  largest  and  distinct  eigenvalues  A| ,  A2, . . . , \k,  and  the  mean  of  Ik*\,  ■  ■  ■ ,  Im  is 
the  sample  estimate  of  the  repeating  lowest  eigenvalue  cr^  of  R.  Following  (11),  the  sample  direction 
estimator  can  be  written  as 

M  K 

Q{<o)=  I  |e:a(a2)P=l-  I  k*Ia(w)|^  (13) 

k=K  +  l  (c  =  l 

The  asymptotic  distribution  of  the  estimates  of  the  eigenvalues  and  eigenvectors  associated  with  the 
distinct  eigenvalues  of  R  is  derived  in  (A.35)-(A.36)  of  Appendix  A.  Corresponding  results  for  the 
conventional  scheme  in  (7a)  can  be  readily  evaluated  as  special  cases  of  this  general  result.  It  is  also 
shown  that  that  the  estimated  eigenvalues  and  a  specific  set  of  corresponding  unnormalized  eigenvectors 
are  asymptotically  (in  the  sense  of  large  N)  jointly  Gaussian  with  means  and  covariances  as  derived  there 
(see  (A.37)  and  (A.38)).  Further,  after  proper  renormalization  and  using  an  exact  relationship  developed 
in  Appendix  B  among  the  different  sets  of  eigenvectors,  it  is  shown  in  Appendix  A  that  (see  A.49) 


£[(?(w)]  =  <?(a2)  +  -  X 


I  .. 

)t  =  i  (A,  -  Afcl 


*  iiki 

/=.  (A.-A,)(A.-A,) 

k*i  1,‘i 


a\w)0J]a(to)  +o(l/N=), 


where  from  (A.18) 

4/,  ^  ^[plR'0kP:R'0j  +  0:r^pJ:r*’0,  +  p:R%r-R^Pk  +  PlR^riyjR'Pkl 


‘  ^2*  '  • »  j  represents  an  M  x  M  diagonal  matrix  with  diagonal  entries  /,,  />,  -  ,  /v.  respectively. 
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In  particular  for  the  conventional  MUSIC  case  with  (A.45)  in  (14),  after  some  simplifications  it  reduces 
to 


£[(?(o;)]=<?(cu)  +  -  I 

N  ,  =  , 


k  =  i (A,  —  At) 

k^i 


+  o(l/N~) 


=  (?((u)  +  ^  i  -  '^'^-,..[(M-K)|/?:a(t^)|--Q(cu)]  +  o(l/N-). 
/V  (A, -<T  r 


(16) 


where  A,,  ;8,,  i  -  1,  2, . . . ,  M,  are  as  defined  in  (7)  and  Q(oj)  is  given  by  (7a). 

Similarly,  from  (A.52) 

Var(<?(w))  A  a'(w) 

^ly  y  y  y  Re[ ( A,„a '( <o)fi,fiJa(co)  +  ))aHLo)p^p\a(w)] 

N  (A,  -AJ(A,-A,) 

k^i  /  / 


+  o(l/N'), 


(17) 


which  again  for  the  conventional  MUSIC  case  reduces  to 
Var((?(w ))  =  cT'iai) 


N  ,  =  , 


AjAfc 


fct,  (A,-Afc) 

k*i 


;\p]a{co)\-\pla{a>f 


XX 

~  ^  2 1)8 la ( <0 )|~i^;a ( w ) P  l  +  o(l/N') 

/  -=  1  VAj  —  Aj } 


=  ^  -^^|/3:a((u)P  +  o(l/N-’). 

N  ,  =  i(A,-o-0 


(18) 


Since  along  the  actual  arrival  angles,  Q(wfc)  =  0,  fc  =  1, 2, . . . ,  K,  (18)  allows  us  to  conclude  that  within 
the  above  approximation, 

cr’(a>J  =  0,  fe=l,2,..., /C,  (19) 

i.e.,  in  all  multiple-target  situations,  where  the  conventional  MUSIC  scheme  is  applicable,  the  variance 
of  the  sample  estimator  for  (7a)  along  the  true  arrival  angles  is  zero  within  a  first-order  approximation. 
It  is  easily  verified  that  this  result  agrees  with  those  of  Kaveh  et  al.,  for  a  two-source  case  (see  equation 
(30)  in  [8]).  An  algebraic  manipulation  shows  that  their  Var(D(wfc))  =  0,  agreeing  with  (19)  here. 

In  particular,  when  all  signals  are  uncorrelated  with  one  another,  that  leads  to  some  interesting  results. 
In  that  case,  the  backward  covariance  matrix  in  (9)  is  identical  to  the  forward  covariance  matrix  in  (6), 
i.e.,  R^=R^-R  implying  equality  of  all  eigenvalues  and  eigenvectors  for  these  matrices.  In  that  case, 
from  (A.20), 

r,ki,  =  '2^^/.^^+p]Y,r]fik).  (20) 
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Let  T7(w)  and  rjica)  denote  the  bias  in  the  f/b  and  the  conventional  null  spectrum  respectively.  From 
(14)-(18)  and  (20), 

iiiai)  =  Q((jj) 

=  S  7~^AiM-K)\p:a{w)\'~Q(w)]  +  o(\/N-)  =  '.rj(w)  +  o{\/N-)  (21) 

2N  .^1  (A,  -or-)- 

and 


—  Q{oj)  X  ^ |/3;a(<i»)|-  +  o( l/N")  =  cr-(a>). 
N  ,_i(A,-o-)- 


(22 


Notice  that  once  again  (r-tw^)  -0,  fc  =  1, 2, . . . ,  K.  Though  there  is  no  significant  advantage  in  using  R 
instead  of  J?'  when  the  exact  covariances  are  known  (except  when  there  is  a  coherent  source),  this  is  no 
longer  true  when  these  covariances  are  estimated  from  the  data.  Then,  in  an  uncorrelated  source  scene, 
only  half  the  number  of  samples  ar**  -equired  to  maintain  the  same  performance  in  terms  of  bias  compared 
with  the  conventional  scheme  while  maintaining  the  same  variance.  Similar  improvements  in  performance 
can  also  be  deduced  for  a  partially  correlated  source  scene  (0<|py|<  1).  To  see  this,  let 


E[u,inuf(t)] 

(£[|«, (t)P]£[|«, (/)!-])'  -’ 


represent  the  correlation  coefficient  between  the  signals  u,(r)  and  u,(t).  Then  in  the  case  of  the  f/b  scheme, 
the  equivalent  correlation  coefficient  p,j  appearing  in  (lOa)  for  the  same  pair  of  signals  can  be  shown  to 
be  [12] 

p„  =  |p,j|  e  ^  cos((  M  -  1  )aj,,  +  <i>,j ),  wi,  -  (oj,  -  w,  )/2. 

Clearly,  Ip.J  ^  lp,y|;  or  stated  in  words,  the  f/b  scheme  has  in  effect  decorrelated  the  signals  beyond  their 
original  correlation  level.  Since  uncorrelated  signals  have  superior  performance  in  this  new  scheme,  any 
amount  of  decorrelation  will  essentially  lead  to  improved  performance. 

The  general  expressions  for  the  bias  and  variance  in  ( 14)  and  (17)  can  be  used  to  determine  the  required 
sample  size  for  a  certain  performance  level  or  a  useful  resolution  criteria.  Though  the  general  cases  are 
intractable,  an  analysis  is  possible  for  two  uncorrelated  sources.  As  shown  in  the  next  section,  the 
performance  of  the  proposed  scheme  can  be  evaluated  in  terms  of  a  resolution  threshold  for  a  two  source 
scene. 


J.2.  Two-source  case 

In  this  section  we  will  assume  that  the  two  sources  present  in  the  scene  are  uncorrelated  with  each 
other.  In  that  case,  from  (21 )  and  (22)  with  K  =2,  the  f/b  scheme  is  uniformly  superior  to  the  conventional 
one  in  terms  of  the  bias  of  the  estimator  by  a  factor  of  two.  This  conclusion  is  also  supported  by  simulation 
results  presented  in  Fig.  1  with  details  as  indicated  there. 

The  deviation  of  plw,)  and  Tj(<u,)  from  zero — their  nominal  value — suggests  the  loss  in  resolution  for 
the  respective  estimators.  Since  the  estimators  have  zero  variance  along  the  two  arrival  angles  in  both 
cases,  for  a  fixed  number  of  samples  a  threshold  in  terms  of  SNR  exists,  below  which  the  two  nulls 
corresponding  to  the  true  arrival  angles  are  no  longer  identifiable.  This  has  led  to  the  definition  of  the 
resolution  threshold  for  two  closely  spaced  sources  as  that  value  of  SNR  at  which  [8] 

£[(?(^<;,)]=£[(?(w:)]  =  £[(?((a;,  +  w,)/2)]. 


(23) 
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Fig.  1.  Bias  at  one  of  the  arrival  angles  vs.  angular  separation  for  two  equipowered  sources  in  an  uncorrelaled  scene.  A  seven 
element  array  is  used  to  collect  signals  in  all  these  cases.  Input  SNR  is  taken  to  be  10  dB. 
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With  the  midangle  =  ( w,  +  a>2)/2,  we  also  have 


!/3la(wn,)|'  = 


l{S>\{MioJ2))- 
1  +  Si(  Mwa) 


(30) 


and 


\p:a(wj\-^0.  (31) 

Using  these  inner  products,  Kaveh  et  al.  found  the  resolution  threshold  ct',  array  output 

SNR)  for  the  conventional  scheme  to  be  [8] 


Nl  A*  1  \  5(M  -2)  /  jj 


(32) 


The  corresponding  threshold  (j  for  the  proposed  (f/b)  scheme  can  also  be  found  by  using  these  norms. 
In  that  case  from  (23), 


iri0(M-2)f  /  2N 

n[  d"  r  (''^5(M-2)'^  )  ■ 


(33) 


This  asymptotic  analysis  is  also  found  to  be  in  agreement  with  the  results  obtained  by  Monte  Carlo 
simulations.  A  typical  case  study  is  reported  in  Table  1.  Under  the  equality  conditions  in  (23),  the 
probability  of  resolution  was  found  to  be  0.3  in  both  cases  there.  This  in  turn  implies  that  the  above 
analysis  should  give  an  approximate  threshold  in  terms  of  ^  for  0.3  probability  of  resolution.  Comparisons 
are  carried  out  in  Fig.  2  using  (32),  (33)  and  simulation  results  from  Table  1  for  0.3  probability  of 
resolution.  Figure  3  shows  a  similar  comparison  for  another  array  length.  In  all  these  cases  the  close 
agreement  between  the  theory  and  simulation  results  is  clearly  evident. 


Table  I 


Resolution  threshold  and  probability  of  resolution  versus 

angular  separation  for  two  equipowered 

sources  in 

uncorrelated  scenes 

Angles  of  arrival 

Uncorrelated 

Uncorrelated  (f/bl 

e,  0. 

Angular  separation 

2to^ 

SNR(dB) 

Prob. 

SNRidB) 

Prob. 

40.00  4.T00 

0,1090 

16 

0.21 

1.1 

0.24 

17 

0..16 

14 

0.18 

18 

0.48 

1.5 

0..19 

19 

0.70 

16 

0.58 

20 

0.81 

17 

0.61 

5.‘;.00  58.00 

0,1.172 

12 

0.20 

9 

0.21 

1.1 

0.40 

10 

0..17 

14 

0.50 

n 

0.45 

15 

0.59 

12 

0.65 

16 

0.8.1 

1.1 

0.67 

120.00  124.00 

0.1860 

6 

0.1. s 

4 

0.25 

7 

0.24 

5 

0..18 

8 

0.48 

6 

0.50 

9 

0..50 

7 

0.57 

10 

0.79 

8 

0.75 

Number  of  sensors  -= 

10,  number  of  snapshots  =  75,  number  of  simulations  - 

100. 

Vi'l  }'.  \(t  1.  ,lulv 


:68 


H.H.  SA  \  Plllai  Arrav  pnwi'win^  fitr  anjih'  of  urrirul 


20 


15 


THRESHOLD 
SNR  (dB) 


f/b  scheme  using  (3.^) 
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conventional  scheme  using  (  52) 
conventional  case  simulation  with  N  =  75 
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Lig.  2  Resolution  threshold  vs.  angular  separation  for  two  equipowered  sources  in  an  unet'rrcl.itcd  sscne  A  ten  element  .irrav  i' 

used  to  receive  signals  in  both  cases. 
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•••  f/b  scheme  using  (.5.5) 

*  f/b  case  simulation  with  N  =  UK) 
conventional  scheme  using  (52) 

♦  conventional  case  simulation  w  ith  N  =  100 
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Fig.  5  Resolution  threshold  vs.  angular  separation  for  two  equipowered  sources  in  an  uncorrelaled  -.cene  A  seven  element  .irrav 

is  used  to  receive  signals  in  both  c.ises. 


The  above  0.3  probability  of  resolution  can  be  explained  by  re-examining  the  arguments  used  in  deri\  ing 
the  resolution  thresholds  (32)  and  (33).  In  fact,  (23)  has  been  justified  by  observing  that  VuriQiw^))  - 

A  /• 

Var(  Q{  to.) )  =  0.  Although  Vai(  QKta,  -F  wO/2))  is  equally  important  in  that  analysis,  using  ( 18 )  and  ( 22 ) 
it  is  easy  to  see  that  these  variances  along  the  midangle  have  nonzero  values.  This  implies  that  though 
and  satisfy  (23),  in  an  actual  set  of  trials  the  estimated  mean  value  of  (?((ra|  -F  wU/2)  will  almost  always 
be  in  the  interval  (0,  2£[(?((w|  +  ai:)/2)])  and  clearly  the  resolution  of  the  two  nulls  in  ()(w)  is  possible 
only  if  this  mean  estimate  lies  in  the  upper  half  of  the  above  interval.  In  the  special  case  of  a  symmetrical 


'signal  IVruessirig 
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density  function  for  the  mean  value  estimate  along  the  midangle,  this  occurs  with  probability  0.5  and  the 
observed  range  may  be  attributed  to  the  skewed  nature  of  the  actual  probability  density  function. 


,  4.  Conclusions 

An  improved  high-resolution  technique  that  incorporates  the  forward  array  together  with  its  complex 
•  conjugated  backward  form  for  estimating  the  directions  ot'  arrival  of  multiple  signals  is  described  here. 

An  asymptotic  analysis  together  with  expressions  for  mean  and  variance  of  this  proposed  f/b  scheme  is 
presented  first;  and  the  conventional  MUSIC'  scheme  is  then  derived  as  a  special  case  of  this  general 
analysis.  In  particular,  when  all  signals  are  uncorrelated,  the  bias  of  the  null  spectrum  estimator  for  the 
f/b  scheme  is  found  to  be  half  of  that  in  the  conventional  case.  Further,  a  resolution  threshold,  which 
depends  on  the  relative  angular  separation,  number  of  sensors,  number  of  snapshots  and  signal-to-noise 
ratios,  for  two  uncorrelated,  equi-power.'d,  closely  spaced  signals  is  derived  here  for  this  new  scheme  and 
this  is  compared  to  a  similar  result  obtained  i”  the  conventional  case  [8].  From  these  comparisons,  to 
detect  two  arrival  angles  under  identical  conditions  one  needs  approximately  half  the  number  of  snapshots 
with  respect  to  the  conventional  case.  These  conclusions  are  also  seen  to  closely  agree  with  the  results 
obtained  from  Monte  Carlo  simulation  studies.  For  this  scheme,  similar  computations  can  be  carried  out 
in  a  two  coherent  source  scene  in  order  to  obtain  the  corresponding  threshold  [12.  I.^]. 


Appendix  A.  Asymptotic  distribution  of  the  sample  eigenvalues  and  eigenvectors 
corresponding  to  distinct  eigenvalues  of  IR 


With  symbols  as  defined  in  the  text,  S'  representing  the  ML,  estimate  of  the  forward/backward  (f  hi 
covariance  matrix  R,  we  have 


{x'(  n  Kx'l  ft  )>  -h  x^l  ft  Hx*'i  ft )) 


ELE  . 


(A.l  1 
(A.:  I 


where 

/r[Sl  =  B.  B  . p,,].  E-ie,.e, . e\,  j. 

.1  -=  diag[A, ,  Ay . Ak.  ff' - (r'l,  diagl  . U.  U  ■  i  ■  ■  ■  f\i\. 

BB  EE  ^  /^,, 

and  B,  E  satisfies  0,,,  e„  '0,  i  -  1, 2, . . . ,  Af  for  uniqueness.  A'-  is  well  known,  the  eigenvectors  are  not 
unique,  and  let  C  represent  yet  another  set  of  eigenvectors  for  S',  i.e., 

S  CLC  .  (A..A) 


where 

C  -  [f, ,  Cy,  .  .  .  ,  Cm  ],  I  A. 4 1 

CC  (A.5) 

For  reasons  that  will  become  apparent  later,  C  is  made  unique  here  by  requiring  that  all  diagonal  elements 
of  Y  given  by 

Y  B  C,  lA.M 


\.o  r  s.'  i 


Iv  I'Jvo 


’70  t)  H  Kwan.  S  I  hllai  Arrav  pnicessth^  for  angle  ol  arm  al 

are  positive  (  v„  '(),  i  -=  1, 2 . M ).  In  what  follows  we  first  derive  the  asymptotic  distribution  of  the  set 

of  sample  eigenv  ectors  and  eigenvalues  of  S  given  by  ( A3)-(  A.6f  and  use  this  to  analyze  the  performance 
of  the  sample  directions  of  arrival  estimator  (?( oi )  in  ( 13 ).  This  is  made  possible  bv  noticing  that  a  I' hough 

the  estimated  eigenvectors  e, ,  . in  (13)  are  structurally  identical  to  their  counterparts 

/i I ,  /3  , .  .  . ,  ,  and  in  particular  have  t’„  ’(),  nevertheless  as  shown  in  Appendix  B,  they  are  related  to 

f,.  i  -  1,2....,  K,  through  a  phase  factor,  i.e., 

e,-c"''c„  I  \,  2 . K  (A. 7) 

and  hence 

K  S  X 

(^(  (t)  I  -  I  '  ^  |e,  a(  (a  )|  -  1  -  ^  [c  ,<i(  <0  )|‘  —  I  —  ^  y,(  w ),  ( A. 8 ) 

I  I  I  I  <1 


vv  here 

V, ( io  I  ■  ,t',  a(  to  o'.  1  1 

fhus,  the  statistical  properties  of  Qiw)  can  be  completely  specified  by  those  of  c,,  i=  1,2 . K.  and 

toward  this  purpv.se,  let 

T-U  .Vt/;  -  1 1.  (A. 10) 

(•  [y'l-A'.' . Kk . A'ol-N  (A. Ill 

7  2  B  Sii  B  CiC  B  -  YLY  .  (A, 12) 

where  V  is  as  defined  in  (.A. 61  with  v„  •().  /-  1,2 . M.  Further,  let 

I  -sA(7'  -  {)  -  -L  V  [',B  (x'(„)(x'(/i))' +xN«)(x''(n)r)B-.t] 

S  N  „  I 

--L  V  [l(:'(„)(;',„))‘ +c''(/i)(;''iu))  )-.{],  (A. 13) 

\  N  II  1 


with 

c'inl  B  .r'(//)  -  \((),  »  K'Wl  (A.14) 

a  nd 

z'in)  B  .V((),  B  R'^B)  (A.l.M 

It  is  easilv  verified  that  these  random  vectors  preserve  the  circular  Gaussian  property  of  the  original  data 
vectors.  .Again,  from  the  independence  of  observations,  asymptotically  every  entry  in  L  is  a  sum  of  a 
large  number  of  independent  random  variables.  Using  the  multivariate  central  limit  thevirem  [1],  the 
limiting  distribution  of  I  tends  to  be  normal  with  means  and  covariances  given  by 

7|»J  L  V  f|',(r',(„)r;*(n)+z^n)c;'*(n))  -AA.]  =  h.  (A. 16) 

V  A'  „  I 


(Here  onwards  whenever  there  is  no  confusion,  we  will  suppress  the  time  index  n)  since 
/:'  I  i(  r ;  I  r ;  I  *  f  >(  r  ^  *  1 1  4 ;  r  [ ',{ x'  ( x'  f  +  x^  .r'- ) '  1  IP,  -  /j :  T  [  S'  ]4,  -  /3 ;  -  A  . 


(A. 17) 
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and 


E[u,ut,]=^  I  [1{E(2UT2[*/,)+  E(zWj*2'l*z^)+  E{z':2';*z[*z^,)  +  E( 

fy  n  =  \ 

Using  the  results^  for  fourth-order  moments  of  jointly  circular  Gaussian  random  variables  and  after  some 
algebraic  manipulations,  we  have 

E[u,ut,]  -  k[p:R%piR'pj + p:R^MiR''0, 

+ p :  j?'r/y;  R^p,+p:  j?‘'y,y;  r^p,  ]  ^  .  (  a.  1 8  > 

In  obtaining  (A.18),  we  have  made  use  of  (A.17)  and  the  fact  that  for  circular  Gaussian  data 

E[2UJ*]  =  P:E[x^x^']Pj  =  p:E[x\xy]yf  =  O.  (  a.  19) 

Here  y^  is  the  inverted  /3*  vector  with  y  „  =  In  an  uncorrelated  source  scene,  it  is  easy  to  show 

that  plRPj  =  y]Ry,  for  all  i,j.  Using  results  of  Appendix  B,  it  then  follows  that  y,  =  i  =  1, 2, . . . ,  /C 
and  [yK  +  i ,  •  ■  • .  iPm]  =  [/?K*r ,  •  •  • .  Pm]K  where  V  is  an  (M  -  K)x(M  -  K)  unitary  matrix.  In  that  case, 

this  together  with  the  identity  p^aica^)  =  0  (or  i  =  K  +  \,  K  +  2.,. . . ,  M  and  Ic  =  1,  2 . K  simplifies  the 

above  expression  into 


Eikij  =  2^,^iiS,kSji  +  pi  yiy*  4* ) . 

Proceeding  as  in  (A.18)  we  also  have 

E[u,,Uk,]  =  \[plR^fi,plR^p,  +  plR'’p,plR*’fij 

+  plR^iki^R^P,  +  plR^yky;R^p,]^  f;,,,. 
Using  (A.12)  together  with  (A. 13)  and  (A. 10)  we  have 

T  =  A+4=U=  yly"  =y{a+-^f]  y\ 

V  N  \  y  N  / 

which  gives  the  identity 


{A.20) 


(A.21) 


(A.22) 


To  d  r  the  asymptotic  properties  of  the  sample  estimates  corresponding  to  the  distinct  eigenvalues 
X, ,  .k^  of  R,  following  [1,6],  we  partition  the  matrices  A,  U,  F  and  Y  as  follows; 


(A.23) 


'  Let  Ti.rj.Zi.r,  be  jointly  circular  Gaussian  random  variables  with  zero  mean.  Then  (refer  to  [15]  E[r|:?r?Z4]  = 
£(Z|Zf  l£[z4zf  1  +  £[z,zU£'(z4?2']- 

Vol  r.  No  July 


t/.zl. 

0  ■ 

,  Y  = 

'Vu 

V,z‘ 

Lt/„ 

t/zzJ’ 

Lo 

Fi. 

.1^2, 

r22. 
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Here  .1,,  and  K,,  are  of  sizes  K  x  K,  etc.  With  i,A.23)  in  (A. 22)  and  after  some  algebraic 

manipulations  and  retaining  only  those  terms  of  order  less  than  or  equal  to  1/v  /V,  we  have 


■  i,  0 

1 

_J - - - 

'V„ 

.  0  cr'/^,  ^  _ 

vW 

1^22. 

+ 


1 

//V 


A,W],  .i,WV 


A,  0 

0  a'Y,.YL 


F,  0 
0  Y..F.YU 


Lw,,.i,  0 


where 

m/,,  =  v7v( 

H-', ;  =  v/7v  V, , ,  W, ,  =  v/N  y, , , 

and  the  column  vectors  h’,  ,  »v,, . . . ,  are  defined  to  be 

^•■1  r 

f  =  U*.  tc 


[  K'l  ,  H', - -  ]  =  W. 


Similarly, 


YY  =  l„  = 


Ik 

0 

'VV,, 

W’i, 

VV'.V 

1 

.0 

Y22YI:. 

vWl 

0 

0  . 

1 

Thus,  asymptotically  for  sufficiently  large  N,  from  (A.28)  and  (A. 24)  we  have 

0=  w,,+ w:„ 

W,,+  V;:»v;:  =  0, 

{/,,  =  vv,,.i,  +  F,+.i,w;„ 

and 

t'.'i  =  -  (r’ W,, . 

Since  S’O,  this  together  with  (A. 25)  and  (A. 29)  implies 

w„=0,  1  =  1,2 . K  and  =  -h  *,  1,7  =  1,  2, . . . ,  A',  i-^j 

which  when  substituted  into  ( A..11 )-( A.32)  gives 
=  u,„  1  =  1,2 . K, 


+  o(l//V),  (A.24) 

(A.25) 

(A.26) 

(A.27) 

(A.28) 

(A.29) 
(A.30) 
(A.31) 

(A.32) 


+  o(  1/  /V  ). 


(a, -A,) 
u„ 


K\,-(r 

From  (A.ll)  and  (A.6)  we  also  have 


/,7  =  1,2 . AC, 

I  =  AC +  1,  AC +2 . M,  7  =  1,2, 


AC. 


G  =  VN(C-B)  =  v' 


/VB(  V-/)  =  v'NB 

L  T,,  Y,  l„ 


(A.33) 


(A.34) 
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which  gives 

or  g,  =  sfN(c,-p,)  =  Bw,,  i  =  \,2 . K. 

This  together  with  (A.  10)  and  (A.33)  gives 

/,  =  A,  +  ( l/v77)/,  =  A,  +  ( 1//N  )U,.,  j  =  1,  2 . K, 


(A.35) 


c,  =  4,  +  {1/v'N)Bh-,  =4,  +  (l/vW)  I  wj„  i=l,2,...,K. 

1 


(A.361 


Thus,  the  estimators  /’  and  c,,  i  =  1,  2 . K,  corresponding  to  the  distinct  eigenvalues  of  R,  are  asymptoti¬ 

cally  multivariate  Gaussian  random  variables/ vectors  with  mean  values  A,  and  /3,,  /  =  1,  2, . . . ,  fC,  respec¬ 
tively.  Furthermore, 


Cov(  /*,  ^  E[u„Ujj]  =  —  r.„..  /,>  =  1,  2, . . . ,  a: 


(A.37) 


1  M  M  1  M  M  r  ^  . 

Cov(c„c,)  =  -I  I  =  I  ‘ 

N  k-i  i^\  Nj,  =  ,  (Aj-Ak)(Aj-A,) 

k  ^  I  I  *j  t  l^J 


(A.38) 


Notice  that  c,  in  ( A.36)  are  not  normalized  vectors,  and  it  may  be  emphasized  that  in  the  case  of  eigenvectors, 
the  above  asymptotic  joint  Gaussian  property  only  holds  good  for  these  specific  sets  of  unnormalized 
sample  estimators.  However,  from  (A.5)  since  the  eigenvectors  c,,  i  =  \  ,2, . . . ,  K,  appearing  in  (A. 8)  are 
normalized  ones,  to  make  use  of  the  explicit  forms  given  by  (A.36)  there,  we  proceed  to  normalize  these 
vectors.  Starting  from  (A.34),  we  have 


1  M 

||c,|l-  =  CiCi  =  1-t-—  I  |h’j,|’>1. 


(A.39) 


and,  hence,  the  corresponding  normalized  eigenvectors  Cj,  i  =  1,2, . . . ,  K  have  the  form 


c.-lk.ir’f.  =  l+TT  I  |h’ 


1  M  1  w 

L  V  1...  |2  I  Is  V 


(iM  \*1^.  I  M  M 

i  ^  t  /  i 


Using  {A.34)  and  (A.18)  we  have 

1  M 

E[c.  1  =  4.-^  E[M]fi.  +  o(  1/  N^) 


(A.40) 


1  r 

=  '-f  -ff,  +  o(l/N-),  /  =  1,2 . K, 

2N  ,  =  ,  (A,  -  A,) 
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since  from  the  asymptotic  joint  normal  distribution  of  these  zero-mean  random  variables  u„(i  ^  j),  their 
odd-order  moments  are  zero.  Thus,  asymptotically  these  normalized  estimates  for  the  eigenvectors 
of  R  are  unbiased  and  the  exact  bias  expressions  are  given  by  ( A.41 ).  Furthermore,  from 

(A.40) 


A  - 

C,Cj  =  1 


M  M  \  ^  ^ 

I  I  k’kjIM  M] 

k=l  I 

k*i  kr,  / _ 


+  7=1  I  ^kAp]+  I  I  I  y^’k.yi'tPkP] 

'  ^  I  »■  I  k  i  I  N  /  I 


I  /MM  ^  ^  M  \i  ^ 

\  k  y  I  I  1 


+  i  I  n,ivi-„i*M;+ 1  I  Wk.\-^tM'i\+o(\/N-). 


k  -•  I  /  .=  1 

k  i*  I  I 


k-|  tk.\ 

k*  t  f*) 


(A.42) 


Again,  neglecting  terms  of  order  1/ N'  and  proceeding  as  above,  this  expression  reduces  to 

,  .  1  /  xr  M  \  1  M  Xf 

e[c.c:]=p,p:-—i  I  £[|H-„n+ 1  e[\wJ]]p.p:+-  i  i  E[w„wt]M',+o(\/N'- 

Ziy  \  k^\  k^i  I  ly  k-\  I  - \ 

\k<^i  k^J  /  k^  t  I*  j 


^  ^  1  w  r . .  ^  r  n 

P  P] - I  .  -  .  ,+  X  -  ,  p,p] 

-  k  ^ I  k*  I 


1  M  M  r 

k*  i  1*1 


and  proceeding  in  a  similar  manner  and  using  (A.21) 


E[c,cJ]  =  l}.l}J-  — 


}  M  r  M  f' 

—  y  ,-h  y  B,pJ 

2N  (A,-AJ'  ,7,  (A,-Ak)' 

_  k  ^  I  k  ^  j  ^ 


1  M  M  /’ 

^  ^  ri  fwi  i  \PkpJ+o(\/N''). 

N  k^\  /-i(A,— Ajij)(A^  —  A|) 

k  ^  ^  ^ 


(A.44) 


An  easy  verification  shows  that  Cov(c,,  c,)  is  once  again  given  by  (A. 38),  but  nevertheless,  ( A.42)-( A.43) 
will  turn  out  to  be  useful  in  computing  the  asymptotic  bias  and  variance  of  the  sample  direction  estimator 
Q(u))  in  (13). 

The  conventional  MUSIC  scheme  [16]  now  follows  as  a  special  case  of  this  analysis.  In  that  case,  from 
( A.16)-(  A.18)  retaining  only  the  terms  corresponding  to  the  forward  array,  we  have 


f  iki,  -  P ,  RPkP I  RPi  ~  A,A,(),/(6ji 


(A.45) 


Signal  Pri>cessing 


B.H.  Kwon,  S.U.  Pittai  /  Array  processing  for  cngle  of  arrival 


275 


Thus,  for  the  conventional  MUSIC  case  using  (A.45),  we  obtain 


1  M  X  X 

E[c.]  =  p.-—  S  — ^A  +  o(l/N^ 

2N  j  =  ,  (A.  -  Aj) 

E[qc,]  =  fifij-—\  I  77— —3+  I  - 


2N 


kt,  (A,-AJ-  (A,-A,r 

i(  I  k^j 


PSl 


1  ^  A, A, 


+  T7  I 


Nkt,  (A,-Ak) 

k  ^  I 


;^,^l6,  +  o(l/N-). 


(A.46) 


(A.47) 


and 


A, A), 


L 


M 

+  I 


AjA)( 


^  =  i(A, 

kiTj 


ppj  -  ^  1  -  5, ) + o(  1/  N=), 


where  Aj,  /3.,  t  =  1, 2, . . . ,  M  are  as  defined  in  (7). 

Once  again  for  the  f/b  case,  using  (A.8)-(A.9)  and  recalling  that  the  eigenvector  estimators  appearing 
there  are  normalized  ones,  we  have 

y,(co)  =  |c;a{w)p  =  aHco)c,c]a{co)  (A.48) 


and  from  (A.43)  and  (13)  we  have 


£[Q(a))]  =  0(cv)  +  -^  I 
N  ,  =  j 


^ iikk 


M 

I  -  - 

aT,  (A, -a,) 
Lk»i 


M  M 

-  I  I  —  -  -  - 

kt,  (A,-A*)(A.-A,) 

i  i 


^ rik; 


7-  a*(co)fi^Ja(w) 


+  o(l/N-), 


(A.49) 


» 


and  this  shows  that 

E[(P((o)]-^  Q((o)  as  N-^oo. 

Also, 


Var((?(a,))  =  E[(?’(a,)]-(£[(?(w)])==c|^^l- I  ) 


K  K 

=  I  I  (E[v,v,]-E[y,]E[y,]). 

r=l y=I 


Using  {A.48)  and  (A.42),  after  a  series  of  manipulations,  we  have 
C[y,y,]-E[y,]E[y,]  =  ^E[d,((u)(/,(a;)]  +  o(l/N^), 


where 


{A.50) 


dfoj) 


It  I  \ 


(u))PxP',o(w)  +  w%a*  (w)Pfila{o}) 


)■ 
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which  gives 

,  2  "  "  Re[|r„„<i’(<»)ftp:o(»)+n,„a'(i,))ftP>(u)la'(u 

£[v,y,]-f[y,]£[y,]  =  -I^  E - *,-*,( - 

I  i*j 

+  o{l/N'). 

Finally,  with  (A.51)  in  (A.50)  we  get 


(A,  A^)(A^  A;) 


,,  .A.  y.  2  ^  ^  «  Rc[{l\,,a\w)p,pla{w)  +  r,,,,afoj)p,fi;a(w)\a\w)i 

Var((?(w))  =  —  I  I  I  I  - - 

TV  ,  =  1  ,  =  I  k  -  1  / -  I 

A;  ^  I  1^1 

+  o(1/N’)-^0  asN-*oc. 

Thus  (?(w)  is  a  consistent  estimator  in  all  cases. 


Appendix  B.  Equivalence  of  eigenvectors 

Let  W  be  an  M  x  M  Hermitian  matrix  with  distinct  eigenvalues  A,,  A;, ... .  A,,  where  m, 

represent  their  repetitions.  Then  m,  +  m.  +  -  ■  •  +  »?,  =  M;  further,  let /3, 1, . fi 

one  set  of  associated  normalized  eigenvectors.  With 

and 

A  =  diag[A  I  ,A|,...,A|,A2,...,Aj,...,A.,...,A,], 
we  now  have 

r  =  BiAb\. 

Let 

=  [4:i  >  P:2 ,  ■  ■  ■ ,  P:m] 

represent  yet  another  set  of  normalized  eigenvectors  of  R.  Then 
B.B:  =  /„,  R  =  BjAbI, 
and  from  (B.l)  and  (B.2)  we  have 
B,Ab^  =  B.ABl 
or  equivalently 

A  V  =  V  i, 

where 

V  = 

Thus,  V  is  also  unitary  and,  further,  .1  and  V  commute.  Moreover,  from  (B.3)  we  have 
A,c,,  =  i',,A,, 

Sign.«l  Procesving 


)pj:a'(w}] 


(A.51) 

%  A 


(A.52) 


,  m,, . . . , 

I  vr,  represent 


(B.l) 


(B.2) 


(B.3) 


(B.4) 
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which  for  A,  A,  gives  =  0.  This  together  with  the  fact  that  V'  is  unitary  implies  V  is  block  unitary  with 
blocks  of  sizes  m,,  m,, . . . ,  m/.  hence  from  (B.4)  we  have 


B:  =  B,  V  =  B 


V, 


where 


1=1,2 . r. 


(B.5) 


Notice  that  in  the  special  case,  when  all  eigenvalues  of  R  are  distinct,  then  V  is  diagonal  and  unitary 
and  each  diagonal  entry  is  a  phase  factor.  In  that  case 

=  i  =  l,2,...,M.  (B.6) 

In  particular,  different  sets  of  signal  subspace  normalized  eigenvectors  of  5  are  related  in  this  fashion 
(see  (A. 7)). 
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